Imaginary Cubic Perturbation: Numerical and Analytic Study 
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The analytic properties of the ground state resonance energy E(g) of the cubic potential are 
investigated as a function of the complex coupling parameter g. We explicitly show that it is possible 
to analytically continue E(g) by means of a resummed strong coupling expansion, to the second 
sheet of the Riemann surface, and we observe a merging of resonance and antiresonance eigenvalues 
at a critical point along the line arg(p) = 5n/4. In addition, we investigate the convergence of the 
resummed weak-coupling expansion in the strong coupling regime, by means of various modifications 
of order-dependent mappings (ODM), that take special properties of the cubic potential into account. 
The various ODM are adapted to different regimes of the coupling constant. We also determine a 
large number of terms of the strong coupling expansion by resumming the weak-coupling expansion 
using the ODM, demonstrating the interpolation between the two regimes made possible by this 
summation method. 
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I. INTRODUCTION 
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with g real positive, has been a subject of a number of 
studies since it has been conjectured by Bessis and Zinn- 
Justin (1992) to have a real spectrum. Arguably, Eq. (JTJ 
provides the simplest example of a VT symmetric Hamil- 
tonian, that is, of a Hamiltonian invariant under a si- 
multaneous complex conjugation (i.e., time reversal) and 
parity transformation x i— > — x (see Ref. [l[). In partic- 
ular, the reality of its spectrum was proven first in the 
strong coupling limit 0, then more generally in Ref. [H 
(together with the positivity of the eigenvalues). 

Here, we continue our investigations Q regarding its 
numerical and analytic properties, including the analytic 
structure of its resonance energy eigenvalues as a function 
of the complex coupling parameter g, specially in the 
strong coupling regime and near the first level crossing 
singularity. 

We denote by E(g) an energy eigenvalue of the Hamil- 
tonian ([I]). One verifies that its perturbative expansion 
contains only integer powers of g: 
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L=0 



with real coefficients El- The series is divergent for all g, 
but is Borcl summablc 5, 6|, and a steepest descent cal- 
culation of the path integral representation of the corre- 
sponding quantum partition function 0413 yields a large 
order behaviour for the nth energy eigenvalue of the form 
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tt is the instanton action. 
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where A 

More recently, Pade summability was also rigorously 
established in Ref. The imaginary part of the en- 

ergy levels on the cut for g = —\g\+i is positive. 
This proof confirms numerical investigations based on 
the summation of the perturbative expansion by Pade 
approximants [12 1. 

More information on analytic and numerical properties 
was provided by an ODM summation [3, [l4| of the per- 
turbative scries for the ground state energy (see Ref. 
It was based on the mapping jh>A. 
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where p is an order-dependent parameter. 

In particular, the leading term of the strong coupling 
expansion could be determined (and was found in agree- 
ment with results coming direct solutions of Schrodinger 
equation [TEl [Tsj ] ) and the positivity of the imaginary 
part verified numerically for all g < 0. It was shown that 
along the negative real axis, the imaginary part above 
the cut is a simple positive decreasing function inter- 
polating smoothly between a strong coupling power-law 
behaviour and a non-analytic exponential tunneling fac- 
tor obtained from semi-classical instanton calculation for 
weak coupling. Finally, it was argued that the ground 
state energy is analytic up to arg(g) = 5ir/4 in the sec- 
ond Riemann sheet, a property that we discuss in this 
article in more detail. 

By a simple coordinate translation, one can relate the 
Hamiltonian ([TJ to another VT symmetric Hamiltonian, 



where \ — 9 
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This correspondence suggests to sum 



the series by another ODM of the form 

A 



(1 - A) 5 / 4 



(6) 
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which is shown here to be much more efficient in the 
strong coupling regime, while the mapping (QJ (and a 
variation thereof discussed below) is more efficient for 
weak coupling. The method © is used here in order to 
determine a number of coefficients of the large g expan- 
sion (small x expansion). Finally, studying the limit of 
weak x, we also construct a corresponding continued frac- 
tion (a subclass of Pade approximants), which enables us 
to obtain more insight into the analytic structure. 

Combining all methods, we have a precise control on 
the analytic continuation of the ground state energy in 
the first Riemann sheet of the uniformization variable 
X = .9~ 4 / 5 , except near the cut in the x-plane. Based 
on our observation, we conjecture that the ground state 
energy of the Hamiltonian ([5]) is a real analytic function, 
with a cut on the real negative x axis, which corresponds 
to a complex phase of arg(g) = bn/A for the original 
coupling parameter g and thus lies on the second sheet of 
the Riemann surface for the original coupling parameter 
g. This domain is not covered by Ref. The imaginary 
part of the ground state energy on the cut for x = ~\x\ + 
i is positive. The small x expansion of the resonance 
energies converges in a disk. From the point of view of 
the initial Hamiltonian (p}, the singularity closest to the 
origin has the interpretation of a level crossing. 

This can be seen as follows: the cut for negative x 
corresponds to a cut of the energy level on the second 
sheet |17l - [f9j of the Riemann surface ( "Bender- Wu cut" ) . 
It starts from the limit of the circle of convergence of 
the strong coupling expansion, i.e., from a point on the 
boundary of the circle of convergence of the x expansion. 
Two resonance energies are attached to points infinitesi- 
mally displaced above and below the cut. The functions 
describing the real parts of the two resonance energies 
are equal on the cut. They cross as we cross the cut. By 
contrast, the functions describing their imaginary parts 
are equal and opposite on the cut and vary smoothly in 
its immediate vicinity. 

At the start of the cut, the imaginary part of the res- 
onance and antiresonancc energy is zero, but the real 
parts still cross, and this configuration therefore corre- 
sponds to a level crossing, given by the confluence of two 
energy levels on the second sheet of the Riemann surface 
characterizing the resonance energies as a function of g. 
For the Hamiltonian ([5]), it corresponds, in addition, to a 
spontaneous breaking of the VT symmetry (namely, the 
spectrum composed of the resonance and antiresonancc 
energies with their equal and opposite imaginary parts 
corresponds to a manifestly non-VT symmetric time evo- 
lution). Starting from the region where the two eigenval- 
ues merge and are real, we can say that after the singu- 
larity, they become complex conjugate. 

To confirm this picture, that is, analyticity in a cut- 
plane in the x variable, sign of the imaginary part of the 
ground state energy on the cut and level merging singu- 
larity we investigate together the ground state and the 
first excited state, forming symmetric combinations: sum 
and difference squared of the corresponding energies. We 



verify that these functions are regular at the level merg- 
ing singularity, showing that the level merging involves 
the ground state and the first excited state. This gives 
us a totally consistent picture and, moreover, allows us 
to determine with increased precision a number of terms 
of the small x expansion. In turn, from the small x ex ~ 
pansion a continued fraction expansion is derived, whose 
coefficients are all positive and seem to converge to a 
value consistent with a square root singularity, within 
errors. 

As a final remark, we notice that the Hamiltonian JT]) 
can be considered as a toy model for the i (f> 3 quantum 
field theory in two and three space dimensions, whose 
renormalization group functions describe the Lee- Yang 
edge singularity of Ising-like statistical models [2(J [2l| . 

This paper is organized as follows. We first discuss the 
scaling of the coupling constant and the related strong 
coupling Hamiltonians in Sec. [H] General properties of 
ODM summation methods are recalled in Sec. IIIII Op- 
timized methods for the problem at hand are introduced 
in Sec. HVl Numerical results for g positive and negative 
are obtained and discussed in Sec. El Sec. IVIl is devoted 
to the strong coupling expansion and level merging. An- 
alytic properties of the ground state resonance in the 
second Riemann sheet and continued fraction expansion 
are discussed in Sec. I VIII Conclusions are relegated to 
Sec. lVTHl 



II. STRONG COUPLING HAMILTONIANS 

The time-independent Schrodinger equation 

[-\^ I + \x^+ l -^jx^^{x)=E{g)^{x) (7) 

corresponds to the Hamiltonian ([T]). With proper bound- 
ary conditions at infinity, this equation determines the 
eigenvalues E{g) of the operator. 

After rescaling the variable x cis x i — y q -^ 5 x, the 
eigenvalue equation can be rewritten as follows (the su- 
perscript [sc] stands for the strong-coupling limit) 

HW(0^(x) = EW(Oi>(x), (8) 

where 

£ = 9~ 2/5 , 0) 

and the strong coupling Hamiltonian is 

The energies of the original and the strong coupling 
Hamiltonian are related by 

E{g) = fl VB £W(j-VS) = ff V6 E ^(0 . (II) 



3 



In particular, we have in leading order 
E(g) ~ g 1/5 E [sc] (0) 

g-s-oo 

and the strong coupling expansion 



(12) 



(13) 



L=0 



In the uniformizing variable -y/f = 9 _1 ' 5 , the singularities 
of E(g) in the complex plane correspond to level cross- 
ings, whereas the cut of the original resonance energy, 
defined as a function of g, is moved to £ = oo. 

We then shift the coordinate x H> x + i £. The equation 
becomes 



Tdx 1 ' 6 



-r + -£ a; 



^ 3 W) = £ Isc] (0 iK*). 



We conclude that 



^ [sc] (e) = -|e 3 + ^ [qc] (e 2 )- 



(14) 



(15) 



where E^ qc '(x = £ ) is a function of the linear coupling 



x = £ 2 = . 9 - 4/5 



(16) 



and thus an eigenvalue of the VT symmetric Hamilto- 
nian ([5]) corresponding to the Schrodinger equation 



where 



(17) 
(18) 



ff[qc] w = -^ +i U^ + ^ 



A generalization of the considerations leading to Eq. ([14]) 
to the normalization of the cubic potential used in 
Ref. [iU is given in Appendix [C] This generalization 
is somewhat non-trivial because the precise form of the 
strong coupling Hamiltonian depends on the normaliza- 
tion of the initial perturbation. When translated to the 
conventions of Ref. [22| . the coefficient of relative order 
£ 3 becomes 1/108 for all states of the spectrum; this af- 
fords an explanation for the observation made originally 
in Table 2 of Ref. H2]. 

The expansion of the ground state resonance energy 
i?[ qc ](x) in powers of \ reads 



oo 



E^\ X ) = Y J VTX 

L=0 



(19) 



Unlike the perturbative expansion in g, it is convergent 
in a disk. Therefore, E(g) has a convergent large g ex- 
pansion of the form 

1 



E(g) = 9 1/5 -o£ 3 + £ [qcI (x) 



1 

3^ 



,1/5 



J2 E 



[qc] g-iL/5 



(20) 



These simple transformations show that E^ (£) has only 
one term odd in £, which moreover can be read off from 
Eq. ([20)) . As a consequence, for g large, the relevant 
expansion variable is g~ 4 ^ 5 rather than g~ 2 ^ 5 , an ob- 
servation that influences the construction of optimized 
ODM summation methods for our problem. In the uni- 
formizing variable g" 1 ^ 5 , it is sufficient to determine the 
eigenvalues for 



< arj 



:(,-«■) 



< 



(21) 



in order to know them in the whole first Ricmann sheet 
of the x variable (— tt < arg(x) < 7r). Conversely, knowl- 
edge of the resonance eigenvalues for — n < arg(x) < it 
implies a determination for — 57r/4 < arg(g) < 57r/4 
which exceeds the first Riemann sheet. 

Indeed, the Hamiltonian (fT7|) is VT symmetric for all 
real values of the "quadratic coupling" x but corresponds 
to the Hamiltonian (JTJ) only in the case x > 0- It i s 
thus interesting to study its spectrum also in the region 
X < 0, which corresponds to the analytic continuation 
g = \g\ exp(5i7r/4). Equation (|17[) is also related to 
residues of poles of solutions of a linear system associated 
to the first Painleve equation [23| . 



III. GENERAL SUMMATION METHODS 
A. Idea of the method 

We here briefly recall the ODM summation method. 
Details can be found elsewhere [III, [l4| . The method is 
based on some a priori knowledge, or educated guess, of 
the analytic properties of the function that is expanded. 
It applies both to convergent and divergent series, al- 
though it is mainly useful in the latter situation. 

In what follows, we consider a function analytic in a 
sector and mappings g h- > A of the form 



9 = P C(A) 



(22) 



where CW is a rea l analytic function increasing on < 
A < 1, such that C(A) = A + C(A 2 ) and, for A -> 1, 
C(A) oc (1 — A)~" with a > 1. A possible choice is 



A 



9 



(1-A)« 



(23) 



L=0 



The parameter a has to be chosen in accordance with the 
analytic properties of the function E (a rational num- 
ber in our examples) and p is an adjustable parameter, 
whose interpretation is an order-dependent, "artificial, 
local radius of convergence" of the truncated, divergent 
perturbative expansion (which becomes smaller as the 
order of the transformation is increased) . An important 
property that singles out relevant mappings here is the 
following: one chooses mappings such that, for g — > oo 
and thus A — > 1, the quantity <? _1 / a has a regular expan- 
sion in powers of 1 — A. A mo re g eneral discussion of the 
method can be found in Ref. jl3l. H^|. 
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After the mapping (|22|) . 25 is given by a Taylor series 
in A of the form 

oo 

E( 9 (\))=Y,Pl(p)\ L , (24) 

L=0 

where the coefficients Pl(p) are polynomials of degree 
L in p. Since the result is formally independent of the 
parameter p, the parameter can be chosen freely. At p 
fixed, the series in A is still divergent, but it has been 
verified on a number of ex amp les (all Borel summable) , 
and proven in certain cases [24| that, by adjusting p order 
by order, one can devise a convergent algorithm. 

The Kth approximant E^ K \g) is constructed in the 
following way: one truncates the expansion at order K 
and chooses p as to cancel the last term. Since Pk{p) 
has K roots (real or complex), one chooses, in general, 
for p the largest possible root (in modulus) p = pk with 
Pr(pk) = 0, for which P' K {p) is small. This leads to a 
sequence of approximants E(g) m E^ K '(g), where 

K 

eW( 9 ) = J2Pl(pk)\ L (9,Pk)\ Pk(pk) = 0. (25) 

L=0 

The term Pk+i(pk) then gives an order of magnitude of 
the error. 

In the case of convergent series, it is expected that pk 
has a non- vanishing limit for K — > oo. By contrast, for 
divergent series it is expected that px goes to zero for 
large K as 

Pk = O (E K 1/K ) . (26) 

The intuitive idea here is that px corresponds to a 'local' 
radius of convergence. 

Let us conclude this section with a few remarks. 
Above, we have used the selection criterion Pk(pk) = 
for the order-dependent pk which are used to calculate 
the polynomials Pl(pk), for L = 0, . . . , K — 1. Alterna- 
tively, one can choose the largest roots px of the polyno- 
mials P'k{pk) = for which Pk is small. Other mixed 
criteria involving a combination of Pk and P' K can also 
be used. Indeed, one empirically observes that the ap- 
proximant is not very sensitive to the precise value of pk , 
within errors. In the ODM method, the determination 
of the sequence of the p^'s is the most time-consuming 
computational step. Indeed, once the pk are known, for 
each value of g, the calculation reduces to inverting the 
mapping (|22j) and simply summing the Taylor series in A 
to the relevant order. 



B. Convergence analysis 

We here give a heuristic analysis [1, [l3| of the conver- 
gence of the ODM method that shows how the conver- 
gence can be optimized. This will justify the choice of the 



class of zeros of the polynomials Pk or P' K and provide a 
quantitative analysis of the corresponding convergence. 

In order to focus the analysis, we assume that we are 
dealing with functions that have the properties of the 
eigenvalues of the Hamiltonian ([1}. We consider only 
real functions analytic in a cut-plane with a cut along 
the real negative axis and a Cauchy representation of the 
form 

E(g) = 25(0) + i- f~ dg> A(g ° , (27) 
K J -co 9(9 - 9) 

where the subtraction ensures the convergence of the in- 
tegral for g — > — oo. 

For the example we discuss here, one can derive by 
a steepest descent calculation @-[sl an asymptotic be- 
haviour of the form 

A( 5) ^oc_ F ^exp(£) J A>0. (28) 
The function E(g) can be expanded in powers of g, 

oo -, ,.0_ j 

E(g) = Y,E L9 \ E L = - -zJrACff), (29) 

the Cauchy representation of El being valid for L > 0. 

The asymptotic behaviour (f28|) then implies a large 
order behaviour of the form ([3]) , 

E L oc (-A)- L T(L + b+l) oc {-A)- L L b L\. (30) 

L— s-oo k— >oo 

A remark is in order. When g becomes negative in 
Eq. ijTj). the expression iy/g becomes real, and the cor- 
responding resonance and antircsonance eigenvalues are 
attached to the upper and lower side of the integration 
contour which extends from g' = — oo to g' = 0. The dis- 
continuity A(g') is the difference of the imaginary parts 
of the resonance and antiresonance eigenvalues. For the 
potential we study here, it has been proved that A(g') is 
positive [ll|, and this is consistent with Eq. (|2"5)) . In ex- 
pression (|27[) , for g > the integrand is positive through- 
out the integration contour (the denominator contains 
a product of two negative quantities). The eigenvalues 
E(g) of the VT symmetric Hamiltonian ([1]) thus are pos- 
itive since 25(0) = 1/2 is positive. Similarly, the expres- 
sion (|29|) shows that the sign of E L is (-1) L+1 , a result 
consistent with the estimate given in Eq. (|30|). 

We now introduce the order-dependent mapping (|22|) . 

9 = PCW, (31) 

with the understanding that the point A = 1 corresponds 
to the point g — > oo. Because we can alternatively for- 
mulate the dispersion relation in terms of A, the Cauchy 
representation then can be written as 



5 



Im(A) 




K factors out in the equation defining the saddle point. 
Thus, A(g) can be replaced by its asymptotic form ([25)) 
for g — > 0_, except for g(X) of order one (or larger) and 
thus A close to the initial endpoint A = 1 for K — > oo. 
The contribution of the endpoint depends on the analytic 
properties of A(g) but is bounded by a quantity of order 
expidK 1 -^ 01 ) [see Eq. gBJ]. 

In what follows we set 



R 



(36) 



FIG. 1: (Color online.) Image in the complex A plane of the 
negative real g axis under the mapping g — p A/(l — A) 5//2 , for 
p = 1/2. A point infinitesimally displaced above the negative 
real axis is mapped onto a point in the upper complex A plane. 



Im(A) 




Refcl) 



FIG. 2: (Color online.) Same as Fig. [T] but for the mapping 
5 = pA/(l-A) 5 / 4 . 



where T initially is the image of the cut on the real neg- 
ative axis (see Fig. Q] and [2] for examples). The contour 
F can then be deformed if the function has analyticity 
properties beyond the first Riemann sheet. 
We expand 



E(g(X)) = £ P K {p) [\{g)\ K 



K=0 



with 



(33) 



(34) 



For K — > oo, the factor \~ K favours small values of |A|, 
but for too small values of A, the exponential decay of 
A(g(A)) takes over. Thus, Pk(pk) can be evaluated by 
the steepest descent method. The ansatz 



R 

PK ~ T7 , R > , 
K-xx> K 



(35) 



implies that in the limit K — » oo, the saddle point val- 
ues of A are independent of K and that g(X) — > 0. The 
former statement can easily been seen by taking the log- 
arithm of the integrand in Eq. (|34l) and observing that 



since p is the only parameter (and it is independent of 
the normalization of g). The behaviour of Pk(pk) is then 
given by the sum of leading saddle contributions. Each 
saddle point contribution is of order 



{Pk{PK )}saddle point = 0(e °) , <J 



MC(A) 



ln(A) , 
(37) 



where A is determined by the saddle point equation 



The analysis is simple only if A at the saddle point is 
real negative. By contrast, if A is complex, the analytic 
properties of A(g) and the possibility of deforming the 
contour T become relevant. 

For p large, the equation has a unique solution (A ^ 
— 1/p), which is real negative. Then we note that at this 
saddle point, as a function of p, 



da 
dp 



1 



C(A) 



>0, 



(39) 



as long as the saddle point value of A is negative. This 
suggests decreasing p as much as possible to improve the 
convergence. The exponential rate a corresponding to 
the saddle point vanishes for 



= Re (ct(A)) = 



1 



mC(a) 



- In(-A) ; 



(40) 



and this defines a special value p = p c of the parameter p. 
In IMEl, we have discussed in some detail mappings 
of the form 



A 



(1 - \) a ' 



(41) 



for various values of a. Most notably, a = 5/2 has been 
applied to the present problem (for a motivation regard- 
ing this mapping, see Appendix |A"| . 

The system of Eqs. (|3"5)) and (|40[) can be solved for the 
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variables p = p c and A = A c . Wc obtain 

• t\ A 
mapping (a): g = p 



(I-A) 5 / 4 ' 
=S> p^ = 3.811 522 ... , \ { c a) = -0.259 901 

• ^ A(l-A/2) 
mapping (b): g = p 



(42a) 



(I-A) 5 / 2 ' 
p(J>) = 4.445 762 ... , A<. 6) = -0.216 262 ... , 

A 



(42b) 



mapping (c): g = p 



(1 - A) 5 / 2 ' 

pi c) = 4.895 690 ... , A< c) = -0.189 645 . 



(42c) 



Then, we expect the error ek at order K to be of the 
same order as the leading saddle point contribution, that 
is, of order Pk(pk) A(g, Pk) K or in view of Eq. (|571) . 

s K = 0(e K ° X(g,p K ) K ). (43) 

For K — > 00, at g fixed, A goes to 1 as 

l-\(g,p K ) ~ CziKg)- 1 ^, (44) 

K— foo 



where for p = p c and methods (|42a[) . (|42b|) (|42c[) . respec- 
tively 



C { 2 a) = {ApP)' 11 * « 10.23, 



(a)x4/5 



c. 



(6) _ / 



(6) 



a 



(c) 



(c) 



2/5 



2/5 



3.40. 



3.53. 



(45a) 
(45b) 

(45c) 



Then, 



A K (3,Pk) ~ exp (-ChK^-^g-yA . (46) 

K — >oo \ / 

Then, two cases are possible: 

Case (i) : if the contribution to the integral correspond- 
ing to <?(A) of order unity decreases exponentially with K 
(which implies the possibility of deforming the contour 
r), we can also choose p < p c and the ODM method con- 
verges for all g ^ 0. Moreover, this in turn implies that 
the function E(g) is an entire function of g~ x l a . This 
clearly is a non-generic situation, but several examples 
have been met in the form of simple integrals. 

Case (ii) : in a generic situation where the mapping re- 
moves the singularity at infinity but other singularities 
are present, as explained in Rcfs. [1, [ID, [l4| we then ex- 
pect the optimal mapping to correspond to p = p c and 
Pk(pk) (or Pk+i(pk) if we choose for px a zero of Pk) 
to be, asymptotically for K —> 00, itself of order 



P K (p K ) = 0[exp (C^K 1 - 1 /* 



(47) 



where a = 5/4 for method (|42a[) and a = 5/2 for meth- 
ods (|42b|) and (|42cj) . respectively. 

The domain of convergence then depends on the sign of 
the constant C3. For C3 > 0, the domain of convergence 
is 



9 < \7T 



cos 



arg(^) 



(48) 



For a = 5/2, this domain extends beyond the first Rie- 
mann sheet and requires analyticity of the function E(g) 
in the corresponding domain. 

For C3 < 0, the domain of convergence is the union of 
the sector | arg(g)| < ira/2 and the domain 



Iff > 



C, 



arg(s) 



(49) 



Again for a = 5/2, this domain extends beyond the first 
Riemann sheet. We will see that this is the situation 
realized in the three methods (|42aj) , (|42b[) and (|42cj) . 

One might now wonder how this analysis is related to 
the zeros of Pjc or Pj^- Choosing the relevant zeros is a 
natural way of finding the region of minimal values. In 
fact, the relevant zeros lie in the region where two leading 
saddle points cancel or the leading saddle point(s) cancel 
the contribution coming from the remaining part of the 
Cauchy integral. 

Indeed, in the discussed examples of case (i), one could 
show that the zeros correspond to a cancellation between 
saddle points. In the case (ii), empirically one observes 
that the zeros are consistent with the asymptotic be- 
haviour (|35p but, in general, with power-law corrections: 



Pk = ^7 (l + S K ), 



(50) 



where 5^ — >• for K — > 00. A comparison with Eq. (f3"T)) 
then leads to the following correction factor to the saddle 
point contribution 



exp 



K S 



K 



Pc C(A C ) 



= exp(K6 K ln(-\ c )) . 



(51) 



Consistency with the preceding analysis, cancellation of 
the endpoint contribution in the integral and convergence 
analysis, then requires 

5 K = OiK" 1 ' 01 ) 



(52) 



and this is also what is generally found. 



IV. OPTIMIZED SUMMATION METHODS 



A. ODM summation with a = 5/4 

We now consider the spectrum of the Hamiltonian ([1} , 
which can be obtained by solving the time-independent 
Schrodinger equation 

~\& + V + x3 ) m = mx) (53) 
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with appropriate boundary conditions. 

The spectrum has a perturbative expansion in inte- 
ger powers of g with real (rational) coefficients alternat- 
ing in sign, of the form ^ , which can be derived up to 
large orders by standard techniques. Instanton calculus, 
based on a steepest descent evaluation of the correspond- 
ing path integral, allows to calculate the classical action 
A relevant for the large order behaviour @. One finds 
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(54) 



According to Eq. (|15j). the strong coupling expansion 
of an cigenenergy E{g) of the cubic contains a coefficient 
gi/5 (-1/3) g- 6 / 5 = -i g- 1 . We can subtract this term 
and introduce the function 



F(g) 



9 E{g) , 



(55) 



which has a regular small g expansion and a large g ex- 
pansion of the form 



F(g)=g 6 / 5 E^(g-±/ 5 )=g 6 / 5 



£ E K 9- iK,b 

K=0 



(56) 

where the coefficients have been defined in Eq. (fl9l) . 
Our aim is to investigate to which extent the strong- 
coupling expansion is described by an ODM summation 
of weak-coupling perturbation theory. To this end, ac- 
cording to the discussion made in Appendix A of Ref. [H , 
we now introduce the mapping (|42a[) which reads 



9 = P 



(1 - A) 5 / 4 



and set 



F(g(X)) = (1 - A)" 3 / 2 cf>(X,p). 



(57) 



(58) 



The function <fi has a regular expansion both at A = 
(g = 0) and A = 1 (g —> oo). The analysis of Sec. |nl] 
then yields p { c a) = 3.811 522 .. . an d thus R = Ap,i a) = 
18.295 306 .. . in the notation of Sec. ImBl [see Eq. (|32"aj) ] . 

Although the summation methods apply to all eigen- 
values, from now on we consider only the ground state 
energy. The first terms of its perturbative expansion are 

1 1 1 930 

E <-" = 2 + 288 "- W^*''' < 59 > 

The convergence analysis of Ref. Q , partially reproduced 
in Sec. IIII1 clearly indicates that in this problem, one 
expects a convergence of the ODM method only in e~ CK<i 
with Q < 1, and this is confirmed here. We find strong 
evidence for a finite distance singularity. It is consistent 
with this property that we can fit the zeros pk of P' K {p) 
using the functional form 



Pk = 



Ap y c 



(«) 



K 



1 - 



12.94 



(A' + 3) 4 / 5 + 11.97 



(60) 



The zeros of Pk can be fitted by an analogous formula but 
with slightly different coefficients for the correction term. 
It has proven convenient to use the latter expression in 
the ODM summation in the place of the exact values. 

There is strong numerical evidence for a singularity 
at g w 0.687 cxp(5i7r/4). (sec Sec. IVUBj) . One thus 
expects the ODM approximants of order K to converge 
to the energy eigenvalues as (see Eq. (|35a|) 



constant x cxp 
consistent with a domain of convergence 



13.8- 10.23 Re (g~ 4/5 ) K 1/5 } , 

(61) 



Re (V 4/5 ) = Rex > -1.351. 



(62) 

Indeed, this expression reproduces well the region where 
the convergence factor becomes close to unity. The do- 
main of convergence contains a segment on the real neg- 
ative x ax is n °t contained in the convergence domain 
obtained by the ODM method (|l2"c|) of Ref. Q , which is 



Re (g~ 2/5 ) =RcVx> -0.0826. 



(63) 



The expression (fSlj) combined with a fit of the relevant 
zeros of the polynomials P' K yields a coefficient of about 
17.4 instead of 13.8 in Eq. (|61[) . Finally, a rough numeri- 
cal analysis of the convergence toward £^ qc ' rather yields 
cxp( — 18 A 1 / 5 ). This implies that either the asymptotic 
regime is reached very slowly, and this is the most likely 
explanation, or our heuristic convergence analysis in this 
particular example is oversimplified. This problem re- 
quires a further analysis. 

Compared with the ODM method of Ref. 0], the con- 
vergence is expected here to be slower at small g and, 
eventually, at extremely large K. At order 150, for g real 
and positive, the transition occurs for \g\ slightly larger 
than 1. Of course, the method with a — 5/2 also con- 
verges in a larger sector, up to |arg(g)| = 57r/4 instead of 
|arg(g)| = 57r/8 for a = 5/4. 

Finally, because the convergence is smooth (unlike 
what happens with a = 5/2 methods) the precision of 
the results can still be improved (see Tables HI [TTJ) by 
using repeatedly the algorithm that to any sequence S n 
substitutes 

Sn i — y / A+2 ~ g "+ 1 . (64) 



S n + SV1+2 — 2S' n +i 

Empirically, we find that it is most advantageous to apply 
this method to the odd and even orders separately. This 
is a variant of Aitken's A 2 process as it is known in the 
literature. Then at order 150, the method (|42a|) remains 
the most efficient one for large coupling \g\ > 1. 



B. ODM summation with a = 5/2 

In Ref. [H , the ODM method with the mapping (|42c|) , 

A 



9 = 



(1 - A) 5 /2 ' 



(65) 
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was applied directly to the energy eigenvalue E(g), which 
has a large-order expansion of the form 



E{g) = g 



1/5 



K=Q 



-2k/ 5 



(66) 



This expansion effectively becomes a large-order expan- 
sion in powers of g~ 4 ^ 5 if one subtracts the term of order 
— 1<7 -6//5 . The latter is the only term not accounted for 

by the powers of g~ 4 ^ 5 - 

The numerical results in Sec. |V] justify the introduc- 
tion of the function F(g), which is defined in Eq. (|56l) . 
and the mapping (|42a[) that precisely takes the special 
form of the large g expansion into account. However, 
while the method (|42aj) is more efficient for g large, the 
method (|42c[) yields more precise results for K very large 
and converges in a larger sector of the complex plane. 

It is thus natural to investigate an ODM method, still 
with a = 5/2, that incorporates in a better way the prop- 
erties of the large g expansion than the method given by 
Eq. (|42c[) . To this end, we again introduce the function 
(|55p and modify the form of the mapping to transfer the 
property of even powers of g~ 2 ^ 5 into even powers of 1 — A. 
A possible choice for a modified mapping is the mapping 
(020), 



-2/5 



A(l- A/2) 
(1 - A) 5 / 2 



(1 - A) [1 - (1 - A) 2 



(67) 



-2/5 



Also, 



2/7 



6/5 



(1-A)- 3 =(JJ [l-(l-A) 

With this mapping, the function 

0(A) = (1-A) 3 F( 3 (A)) 



21-6/5 



(68) 



(69) 



is indeed only a function of (1 — A) 2 . However, in contrast 
to the previous mapping, the symmetry g 



-2/5 



-fJ 



-2/5 



is only implemented approximativcly, at a finite order, 
because the series in powers of A is truncated. 

Since one cannot expect geometric convergence of the 
ODM method defined by Eq. (|42b>|) . the relevant value 
of \i is n { c b) = 4.445 762 607 .. . (see Sec. IhTbI) . It follows 



that C?> = (i?/2) 2 / 5 = (^ b 72) 2 / 5 = 3.401 535 .. . [see 
Eq. (pB)]. 

For the ground state energy, we empirically find that 
the location of the zeros of P' K can be fitted by 



PK 



A \x c 



1 



5.0 



(K + 2) 4 /s + 4.6 



(70) 



One also expects a term decreasing only like K~ 2 / 5 
but its coefficient is apparently too small to be de- 
tected compared to the larger K~ 4 / 5 contribution. The 



same remark actually applies to the method (|42cl) . We 
have not compared the convergence of methods (|42b|) to 
method (|42c|) systematically, but methods (|42b[) is defini- 
tively better for g — > oo than with the initial a = 5/2 
mapping of Eq. (|42c[) (see Tables HI and ITT|) . For example, 
at order 55, with the new mapping (|42b[) . one finds the 

value E^ = 0.372 545 (9) for the leading strong-coupling 
coefficient defined according to Eq. (|T9")) , corresponding 
to a relative error of the order of 10 -7 , that is, two order 
of magnitudes better as compared to the mapping (|42c[) . 
At order 150, one finally obtains E$ c] = 0.372 545 78(6), 
a value that differs by 4 x 10~ 9 from the strong coupling 
limit in ([9"2"|l . 

Values of the ground state energy E(g) for various real, 
positive and negative values of g, and order 55 and 150, 
are reported in Tables HI— IIV1 

The general conclusion is: at the order we calculate, 
as expected, the results for the modified method (|42b[) 
are more precise than for the more straightforward 
method (|42c|) if g is large, but comparable or worse for 
g small. However, for \g\ > 1, results given by the 
method (|42a[) with a = 5/4 are always the best. This 
is consistent with the fact that higher-order corrections 
of the form g~ 4K / 5 are well represented by the corre- 
sponding expansion for E^ qc \g). Note that the correc- 
tion of relative order g~ 6 ^ 5 for E^ sc \g) provides for a 
certain irregularity in the expansion (|13|) . As a conclu- 
sion, the advantages of the modified version (|42b[) seem 
to be somewhat limited, and we have not investigated 
the method more thoroughly, but we have included the 
results in the tables to stress the general consistency of 
the different implementations of the ODM summation 
methods. 



NUMERICAL RESULTS: GROUND STATE 
ENERGY 



We first display a few typical results obtained for g 
finite, positive and negative, by the various methods we 
discuss in this article: the ODM method with a = 5/4, 
the methods of Sec. IIVBI and Ref. Q and, finally, the 
continued fraction of Sec. I VII B I below. 



A. Real positive axis 

For g finite and real positive, we have chosen the same 
values as in Ref. Q and added results for the case g = 0.5. 
Most results are displayed in Table U for the results at 
order 55, to allow for a direct comparison with Q and 
give an idea about the rate of convergence, and in Table [TT1 
at order 150. We use everywhere the convention that the 
numerical uncertainty applies to the digit in parentheses 
(this may also affect the preceding digit by a shift ±1). 

For g = 0.5, at order 150, with the ODM method given 
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TABLE I: The ODM mappings given in Eq. (|42a[) . (|42b|l and (|42c|) are compared with respect to their convergence properties 
on the real positive axis at a transformation order of 55. Method (|42a[l with a — 5/4 and extrapolation clearly provides the 
best answers for large and moderate coupling. The modified method (|42b|) with a = 5/2 interpolates between method (I42al) 
and (|42cl) ; indeed, the application of the straightforward ODM method (I42c[) still provides for the best convergence at weak 
coupling. 
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0.5 


1 


Method of Eq. (|4"2a)l 
Method of Eq. (|42b|) 
Method of Eq. (|42c) 


0.51689 17642 53171 97821 11588 95(6) 
0.51689 17642 53171 97821 11588(9) 
0.51689 17642 53171 97821 11588 9(6) 


0.53078 17593 04176 67113 556(1) 
0.53078 17593 04176 67113(4) 
0.53078 17593 04176 671135(6) 
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5.0 


21.6 


Method of Eq. 
Method of Eq. 
Method of Eq. 


(|42b)l 
(|42cTl 


0.60168 39332 05191 96159(1) 
0.60168 39332 05(2) 
0.60168 39332 05(2) 


0.73340 99204 85427 96(4) 
0.73340 9920(5) 
0.73340 992(1) 



TABLE II: Same as Table [T] but for a transformation order of 150. For weak and moderate coupling, the method 
(|42cp now takes over. For strong coupling, the convergence of method (|42a|l with a — 5/4 clearly is superior to 
that of the other ODM methods. For completeness, we give also the results obtained from the simple continued 
fraction method (|99jl . For the method defined by (|42ap . the numerical results are post-accelerated using the 
algorithm (|64[) . 
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0.5 



1.0 



Method of Eq. (|42a)l 
Method of Eq. (|42b|) 
Method of Eq. (|42c|) 
Continued fraction 



0.51689 17642 53171 97821 11588 9 
0.51689 17642 53171 97821 11588 9 
0.51689 17642 53171 97821 11588 9 
0.51689 17642 53171 97821 1(0) 



0.53078 17593 04176 67113 55618 18032 
0.53078 17593 04176 67113 55618 18032 
0.53078 17593 04176 67113 55618 18032 
0.53078 17593 04176 67113 55(7) 
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5.0 



21.6 



Method of Eq. (|42a) 
Method of Eq. (|42b|l 
Method of Eq. (|42c|) 
Continued fraction 



0.60168 39332 05191 96158 93564 94(4) 
0.60168 39332 05191 96159(0) 
0.60168 39332 05191 96158(9) 
0.60168 39332 05191 96158 936(0) 



0.73340 99204 85427 96459 24020(0) 

0.73340 99204 8542(8) 

0.73340 99204 854(3) 

0.73340 99204 85427 96459 240(3) 



in Eq. (|42b[) , one finds 

£(0.5) = 0.51689 17642 53171 97821 11588 

956621775(8), (71) 

and with the method given in Eq. (I42cl) 

£(0.5) = 0.51689 17642 53171 97821 11588 

95662 17760 99999 61207 (4) , (72) 

the error being in agreement with the estimate given in 
Ref. Q. Finally, with the method given Eq. (|42a[) one 
finds 

£7(0.5) = 0.51689 17642 53171 97821 11588 

956621775(8) (73) 



and after extrapolation using the method ([M]). 

£(0.5) = 0.51689 17642 53171 97821 11588 

9566217760 99999 612(1). (74) 

For <7 = 1, from a numerical solution of the Schrodingcr 
equation one obtains 

£(1.0) = 0.5307817593 04176 67113 

5561818032 225. (75) 

At order 150, with the mapping (|42c[) . the result reads 

£(1.0) = 0.5307817593 04176 67113 

5561818032 22595(1). (76) 
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At order 150, the method (|42b|) yields a result with 
an error of 10 -33 , which is slightly less precise. The 
method (|42a|) together with extrapolation (J64j yields 

£(1.0) = 0.5307817593 04176 67113 5561818032 

225951(1). (77) 

Clearly, for small values of g, that is, g < 1, at order 150 
the method of 0] gives the most precise results. The en- 
tries for E(1.0) exceed the maximum number of columns 
available in Table [TTl for results obtained at transforma- 
tion order 55, see Table HI 

Finally, to compare with the Pade summation of [l2j , 
we give results for g = 288/49, which corresponds to 
A 2 = 1/7 in [lj]. In Ref. [13], at order 192, the reported 
result is E = 0.61273 81063 88986. In Ref. Q, using the 
method (|42c[) . at order 55, one finds 

) = 0.61273 81063 89(1) . (78) 

With a = 5/4, using method (|42a|) and convergence ac- 
celeration, at order 55, one obtains 

£(lf ) = 0.61273 81063 88984 124(7) . (79) 

This improves over the ODM method @23 in Ref. i by 
six orders of magnitude. Then, at order 150, the result 
becomes 

E{^§) = 0.61273 81063 8898412476 20895 52(6). (80) 

With the method (|42cp . at the same order 150, one ob- 
tains 

£(lf ) = 0.61273 81063 8898412476(3) . (81) 
We have also verified that 

E(^§) = 0.61273 81063 8898412476... (82) 

by a numerical solution of the Schrodinger equation. 

The general conclusion is: for values g < 1, at order 
150 the method with the mapping (|42c[) of Ref. [H is the 
most precise, and the advantage increases with increasing 
order and decreasing parameter g. By contrast, for g > 1, 
at least up to order 150, the method with a — 5/4 takes 
over. Note that from now on and, in particular in the 
tables, we only quote the results obtained by the method 
(|42ap after extrapolation. 

Finally, the continued fraction of Sec. IVII B[ con- 
structed from the strong coupling expansion (but which 
incorporates additional information about level merg- 
ing), provides a rather good representation of the func- 
tion in a wide domain. 

For completeness, in the spirit of reference [l2|, we give 
some indication about the summation of the perturbativc 
series by a continued fraction expansion. We define (for 
the ground state) 

h o (g)=E(g)/E(0), (83) 



and introduce the relation 

Vi(.9) = l + ^, h p (0) = l, (84) 

which allows calculating h p from and determining 

recursively the coefficients k p . The truncated contin- 
ued fraction (obtained by replacing, at some order p, 
hp(g) by 1) generates, alternatively, [n + 1/n] and [n/n] 
Pade approximants. For a Stieltjes function, all coef- 
ficients k p are positive. Moreover, for a divergent se- 
ries with a large order behaviour of the form Q one 
expects the coefficients n p to grow asymptotically lin- 
early with p and the error for the continu ed f raction 
truncated at order p to behave like exp(— C\Jpjg). This 
is indeed what is observed. The coefficient of the lin- 
ear term is compatible with 5/48. More precisely, a 
good fit is k p = (lOp + 3 x (-1)p)/96 + 0(1/ y/p) for 

p — > oo. Also, the coefficient of \fp~fg~ extrapolates to 
C = 2-^/48/5 = 6.19... with a good precision. Com- 
pared to the various ODM summations, the convergence 
is limited to the first Riemann sheet, is poorer than the 
ODM method with a — 5/2 for g small and is much 
poorer than the method with a = 5/4 for g large. For 
example, at order 150, for g = 0.5 the error is about 
1.5 x 10~ 42 , for g = 1 about 6 x 10 -3X , for g = 5 about 
10~ 14 and about 2 x 10~ 7 for g = 21.6. 

B. Negative real axis 

Results for four values of g < are displayed in Ta- 
bles IIIII and IIVI From the numerical evidence, we con- 
clude that ODM methods with a = 5/2 converge also 
for g < but, from the theoretical analysis, the con- 
vergence is expected to be poorer than for g > 0, in 
agreement with the data. The ODM method (|42b[) with 
a = 5/4 converges also on the real negative axis for \g\ 
large enough, that is, \g\ > 0.53 approximately. 

From the analysis of Ref. [|| , we know that the imagi- 
nary part itself is an analytic function with singularities 
on the real negative axis only at g = and at infin- 
ity. Moreover, it is a simple positive decreasing function. 
Compared with the method with a = 5/2, the positivity 
of the imaginary part for g = — \g\ + iO can again be ver- 
ified, but with higher precision for g large. For \g\ — > oo 
and Im (<?) = 0+ , from the leading term in Eq. (p?2"j) suit- 
ably rotated into the complex plane, one now infers 

E(g)/\g\ 1/5 = 0.3013958756586835717823(7) 

+ 0.2189769214314493762936(0) i , (85) 

a result that uses the perturbative series up order 150. 

For g = -0.5, the method (|42b[) with a = 5/4 is not 
expected to converge and, indeed, numerical indications 
are that this is the case. The result is somewhat equiv- 
alent to a direct summation of the initial asymptotic se- 
ries, with an error of the order of the imaginary part 
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TABLE III: On the real negative axis, at order 150, the real part Re E(g) is approximated by several variants of ODM 
methods. The straightforward ODM method (|42c[) provides the best results for weak coupling, whereas method (|42a| 
with a — 5/4, post-accelerated, is superior in the moderate and strong coupling domain. Note that the simple continued 
fraction method (|99l) also provides for a well converged answer. 



0.5 



1.0 



Method of Eq. (|42aj) 
Method of Eq. ([42b]) 
Method of Eq. (file)) 
Continued fraction 



0.47642 7408(3) 
0.47642 74083 2(7) 
0.47642 74083 27179(5) 
0.47642 74083 271(9) 



0.44252 0045124688(4) 
0.44252 004512(4) 
0.44252 00451 24(7) 
0.44252 00451 24688 3662(3) 



5.0 



21.6 



0.43389 06678 10363 12813 1(1) 0.55405 35184 61013 80317 898(0) 

0.43389 06678(1) 0.55405 3518(5) 

0.43389 0667(9) 0.55405 35(2) 

0.43389 06678 10363 12813 116(9) 0.55405 35184 61013 80317 898(0) 



Method of Eq. (|42aj) 
Method of Eq. (j42b|) 
Method of Eq. (p2c| 
Continued fraction 



TABLE IV: The imaginary part Im E(g) at transformation order 150 shows the same characteristics as the real part 
discussed in Table [TTT1 



0.5 



1.0 



Method of Eq. (|32aj) 
Method of Eq. (|42bj) 
Method of Eq. (pEc)) 
Continued fraction 



0.00026 6662(1) 
0.00026 66618 8(2) 
0.00026 66618 82408(1) 
0.00026 66618 824(6) 



0.0155179258 2059(4) 
0.0155179258 2(0) 
0.0155179258 20(6) 
0.01551 79258 20594 2572(2) 



5.0 



21.6 



Method of Eq. (|42aj) 
Method of Eq. (j42b|) 
Method of Eq. (p2c|) 
Continued fraction 



0.18385 80861 86171 172893(3) 
0.18385 80861(9) 
0.18385 8086(0) 

0.18385 80861 86171 17289 33(1) 



0.35140 17775 93691 93624 451(6) 
0.35140 1777(6) 
0.35140 18(0) 

0.35140 17775 93691 93624 45(1) 



on the cut. However, compared to the direct summa- 
tion of the series, in some range the successive results 
first oscillate around the exact value with the order K 
with a rather slowing increasing amplitude rather than 
blowing up. Otherwise, we notice a situation comparable 
to the real positive axis. For \g\ < 1, the ODM method 
of [H yields the most precise results, while for \g\ > 1, the 
method (|42aj) with a = 5/4 and convergence acceleration 
takes over. However, the continued fraction of Sec. lVIIBl 
gives the most precise results in the intermediate range 
and is equivalent for \g\ large. 



VI. STRONG COUPLING EXPANSION AND 
LEVEL CROSSING 

We here demonstrate that considerable information on 
sublcading corrections to the strong coupling asymptotics 
can be obtained by investigating the ODM resummcd 
weak-coupling expansion of the energy levels. 

In order to explore more thoroughly the strong cou- 
pling expansion and the analytic properties of the ground 
state energy of the cubic Hamiltonian, it is very useful to 
also consider the linear coupling Hamiltonian <JTTJ) , 

fiM =-^ +i 3+ H- (86) 
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because a perturbative expansion in the variable 

x = g- A/ \ (87) 

of the energy level of the Hamiltonian H^ qc \ defined in 
Eq. (1171) . is equivalent to a large g expansion of the corre- 
sponding energy of the Hamiltonian ([1]) . We recall here 
that the two Hamiltonians are only equivalent for \ > 0, 
the negative x axis corresponding to arg((/) ~ 57r/4. 

For real g, the ground state energy E(g) of the Hamil- 
tonian ((T|) is a real analytic function. From the numerical 
evidence provided by the ODM summation in Ref . [|| , we 
conjecture that, that, in the variable g -1 / 5 , the eigenval- 
ues can be continued up to | arg (g -1 / 5 ) | = 7r/4. This 

implies that in the variable x> the eigenvalues E^ qc \x) 
are real analytic with singularities only on the negative 
real axis. The series in powers of x are convergent in 
a disk. In the case of the ground state, at the singu- 
larity nearest to the origin (numerical results indicate 
X = Xc = —1.3510...). Note that, by contrast with a 
Hcrmitian Hamiltonian, VT symmetric Hamiltonians do 
not experience eigenvalue repulsion, which explains why 
it is possible for eigenvalues to merge for x real. 



A. Strong coupling expansion 

We now concentrate on the behaviour of E(g) for g — > 
oo, or, equivalently, on the small % expansion given in 
Eq. (THD of £ [qcl (x): 

fi H W = E^/' (88) 

K 

First, we have calculated the values of the coefficients of 
the large coupling expansion as determined by the ODM 
method (|42a|) with a = 5/4, with the improvement by the 
algorithm (fM)) . The successive coefficients of the small x 
expansion are related to 4>(X, p) and its derivatives taken 
at A = 1 (sec Appendix lAl for the first terms). At leading 
order, one finds 

0(1,^=^/5^0]. (8Q) 

For the leading term, we obtain 

E [ ^ c] = 0.37254579045220709825 0601(1). 

Note that even at order 150, with the ODM method (|4"2"c|l 
for the coefficient i? qc ' one still finds an uncertainty of 

4.0 x 10~ 6 . With the method (|42b) . one obtains E^ c] = 
0.372 545 78(9). 

From the ODM method l|4"2"a|) with a = 5/4 applied 
to the perturbative expansion of the function (|5o| . at 
order 150, we have then determined, using the algo- 
rithm described in Appendix [Bl with decreasing rela- 
tive precision (about 10~ 6 for the last term), all terms 
up to order K — 20. The analysis of the behaviour 



of the coefficients with increasing order strongly sug- 
gests the existence of a square root singularity located 
at x = Xc = —1-351 ± 0.002. A square root singularity 
is consistent with the existence of a level merging. To 
confirm this analytic structure and locate the singularity 
more precisely, we have assumed that the point Xc cor- 
responds to a merging between the ground state energy 
and the first excited state energy. 

B. Level merging and strong coupling expansion 

In addition to the ground state, as described, we have 
also generated the perturbative series for the first ex- 
cited state. We do not report here the details of the 
numerical study of the energy of the first excited state. 
The behaviour of the series with respect to the ODM 
summation methods is very similar. The coefficients of 
the strong coupling expansion can be determined with 
comparable precision. The convergence properties of the 
ODM method (|42cp indicate that, again, the eigenvalue 
is analytic in a cut-plane and that the first singularity 
corresponds to the level merging with the ground state 
energy. 

We have then formed the two symmetric combinations, 
the half sum Sqi and the half difference squared Aoi (GS 
and ES stand for ground state and first excited state, 
respectively) : 

A — 1 ( p[ qc l ph c A 2 <; — 1 ( ipM i pI^A 

A 1 - 4 (^ Gs - E ES ) , boi - ^ ^GS + ^ES ) ■ 

(90) 

If the ansatz is correct, the two functions are not singular 
at Xc and Aoi must vanish linearly at Xc- This is indeed 
what is observed (see Figs. [3] and H]). The result confirms 
that at Xc the eigenvalues corresponding to the ground 
state and the first excited state, merge and for x < Xc 
become complex conjugate. 

More precisely, the direct summation of the series for 
A i with the ODM method of Sec. IIVAI shows that A i 
vanishes linearly with x at the point 

Xc = -1.3510415966(3) , (91) 

a result fully consistent with the direct study of 
the ground state. At this point, E^ qc \xc) = 
0.41330579447(3). As a by-product, one also obtains 
for the ground state, the value of E^(x = — 2 4 / 5 ) = 
0.38985020(5) - 0.364427(9)i, which confirms the result 
coming from the continued fraction in Table IVII 

By the same method as for the ground state, we have 
determined the strong coupling expansions of Sqi and 
Aoi- From these strong coupling expansions, one can 
recalculate the strong coupling expansion of the ground 
state energy. The results are completely consistent with 
those from the direct expansion. However, the precision 
is improved for the higher order coefficients as expected 
since the singularity at Xc is now explicitly generated. 

As a necessary ingredient for a more precise deter- 
mination of the coefficients of the continued fraction of 
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Sec. IVIIBI for the ground state, we have calculated the 
strong-coupling expansions of the difference A i and the 
sum Sqi up to to order 28. We have then inferred the 
strong coupling expansion of the ground state at the same 
order although we here give only the first 20 terms. These 
are useful for reference purposes and read as follows: 



E lqc - 
+ 0. 
+ 0. 

- 0. 
+ 0. 

- 0. 
+ 0. 
-0. 
+ 0. 
-0. 
+ 0. 
-0. 
+ 0. 

- 7. 
+ 5. 

- 3. 
+ 2. 

- 1. 
+ 1. 
-0. 
+ 0. 



(x) = 0.37254 57904 52207 09825 06011(5) 
36753 58055 44193 60353 04(6) X 
14378 77004 15066 51583 39(0) X 2 
02658 61056 27059 38713 52(9) x 3 
00988 71650 79200 88729 05(5) X 4 
00461 00192 93623 15160 2(3) x 5 
00240 93426 35048 47521 1 (7) x 6 
00134 88515 29319 85498(8) X 7 
00079 06119 76816 97837(2) x 8 
00047 88478 41414 5725(4) x 9 
00029 72375 58426 7145(5) x 10 
00018 80657 95326 713(9) X n 
00012 08255 49560 587(6) X 12 
86045 58627 946(5) x 10~ 5 x 13 
16744 64642 199(1) x 10~ 5 X u 
42729 47828 030(3) x 10" 
29050 20869 87(5) x 10~ 5 
5409159219 76(3) x 10" 5 X 
04266 03452 042(4) x 10" 
70913 87535 56(4) x 10~ 5 A 
48450 81660 66(0) x 10~ 5 X 20 + 0( X 21 ] 



5 x 15 



17 



5 x 18 

,A9 



(92) 



Note that the errors are strongly correlated. 

This expansion is also consistent with the 20 first coef- 
ficients reported in Ref . fl5| , which have been determined 
with a 10 -10 -10 -11 relative precision by a numerical so- 
lution of the eigenvalue equation. It is also consistent 
with results obtained for the few first terms from a nu- 
merical solution of the Schrodinger equation [l6j |. which 
have relative errors of order 10 ~ 9 . 



VII. CONTINUED FRACTION AND 
ANALYTIC PROPERTIES 

A. General remarks 

To discuss more thoroughly the analytic properties of 
the ground state energy, it is convenient to consider the 
Hamiltonian (fTTj) and parameterize the energy in terms 
of the coupling X - 

Inverting the relations (fTTj) and (|19l) . one finds a form 



relevant for the large X behaviour of E^ (x) : 

i? [qcl (x) = ^ 3/2 + X 1/4 i?(x- 5/4 ). (93) 

The perturbative expansion of E(g) in g translates into a 
large X expansion, because large X corresponds to weak 
coupling g, 

x l/^ (x -5/4 )= l x l/4 + _^ x -l 



155 
13824 



X 



-9/4 



In particular, the imaginary part for X 
the expansion 



(94) 
iO has 



Im #%) 



155 
13824^2 



1 



2V2 M 

9/4 



1/4 



X 



(95) 



For |x| 3> 1, the imaginary part is clearly negative. The 
question that arises is its sign on the cut for |x| of order 
unity. 

Note that the large X behaviour (|M|) implies that the 
Cauchy representation for the perturbative expansion co- 



efficient E [ P can be written as 



El 



[qc] _ 1 [** Im EM( X ) 



X 



K+l 



dx. 



(96) 



converges only for K > 1. The conjecture Im E^( X ) < 
0, which will be substantiated in Sec. IVIIBI is thus con- 
sistent with the signs of coefficients of the expansion (|92p . 
Finally, some additional information can be obtained di- 
rectly from the ODM method with a = 5/2 (see Ta- 
ble IVI[) and also by summing the small X expansion (|92|) 
as we will show. 



B. Continued fraction 

To continue the expansion (p?2"j) outside the circle of 
convergence, we introduce the continued fraction expan- 
sion of 



£ [qc] (x) 



£ [qc] (x) - £ [qc] (0) 
£[qc]'(o) x 



(97) 



motivated by the prejudice that Et qc ' is a Stieltjes func- 
tion since 



Im -B [qc] (x) 



Im^[ qc l(x) 
£ [qcl '(0) X 



(98) 



is positive if Im£g qc ' is negative and i?[ qc l behaves like 
X 1 / 2 for |x| — > oo. We then construct the continued 
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TABLE V: Coefficients a p of the continued fraction for K — 150. The singularity is consistent with a limit 0.1850424 . . .. 



p = l 


p = 2 


P = 3 


p = 4 


0.39122 09320 72635 98993 


0. 18489 83296 22869 56168 


0.18699 38616 55333 76095 


0.18768 40668 9188109149 


p = 5 


p = 6 


P = 7 


P = 8 


0.1851910686 95077 9010(9) 0.18477 4976141944 774(1) 


0.18470 2224817836 84(5) 


0.18507 38286 03338 2(9) 


p = 9 


p = 10 


p = 11 


p = 12 


0.18513 42403 94894(8) 


0.18510 363919575(7) 


0.18499 98937 209(1) 


0.1850158945 34(4) 


p = 13 


p = 14 


p= 15 


p= 16 


0.1850417233 0(4) 


0.18506 43232 9(6) 


0.0.18504 19101(5) 


0.18503 66878(1) 



fraction expansion for the function £?[ qc ] by introducing 
the recursion relation 

fp-iti) = ! + /p(0) = 1, (99) 

mx) 

with the initial condition /o(x) = i^ qc ](x). It allows us 
to calculate the coefficients a p recursively. The truncated 
continued fractions, obtained by replacing at some order 
fpix) by 1, generate, alternatively, [n + l/n] and [n/n] 
Pade approximants, depending on whether p is odd or 
even. 

For a Stieltjes function, all coefficients a p are positive. 
Indeed, using the coefficients ([92]) . which are obtained 
from the series at order 150, we find a p > (see Table W} 
for p < 27 and the coefficients seem to converge slowly, 
within errors, toward the limit expected for a square root 
singularity at x = Xc- V( 4 Xc) = -0.18504241514.... 
For example, the two last coefficients determined with 
sufficient precision are 024 = 0.1850464(8) and a 2 5 = 
0.185042(8). 

To test the idea further, we have substituted in the 
infinite continued fraction 

/taM^ \ + \V^ix, (100) 

Pmax varying from 9 to 27. This amounts to taking a p 
constant, a p — 7 w 0.1850 for p > p max . Remarkably 
enough, for x ~ * 00 1 the corresponding approximation for 
i?[ qc l(x) behaves like x 3 ^ 2 with a coefficient close to the 
exact value 1/3. We have thus adjusted 7 more precisely 
to get 1/3 exactly. This yields a very stable sequence up 
to p = 25, converging slowly toward the expected limit 
0.18504...: 

7i 7 = 0.1850687 , j is = 0.1850261 , 719 = 0.1850603 , 

720 = 0.1850239 , 72i = 0.1850562 , 722 = 0.1850310 , 

723 = 0.1850562 , 724 = 0.1850367 , 725 = 0.185048(4). 

(101) 

corresponding to a singularity located at x — Xc ~ 
— 1.3504, consistent with a direct analysis of the be- 
haviour of the coefficients in the expansion (f9"2")l . As a 



consequence, in all Tables and Figures we have reported 
results obtained from this approximated continued frac- 
tion. 



C. Spectrum for x < 

The convergence of the various methods is rather poor 
for x < 0, specially in the neighbourhood of the cut. 
Table IVT1 displays a few results for i?[ qc l(x). According to 
Table fVTl the most precise results are in general obtained 
from the modified continued fraction, with a precision 
that for I g | > 1 is comparable with the ODM method 
defined by Eq. (|42a[) . Note that for the first value, which 
corresponds to \g\ = 0.5, the ODM method with a = 5/4 
does not converge, as expected, and does not even provide 
a reliable estimate. The method of Scc. lIVBl is not useful 
either. Finally, in the neighbourhood of the singularity 
X = Xc a calculation based on the determination of A01 
and 5*01 gives the most precise results (see Figs. [3] and U]). 

The imaginary part after the cut for x = — |xl + i0 is 
necessarily negative when calculated from the approxi- 
mated continued fraction. Therefore, the general consis- 
tency of the results coming from the continued fraction 
and the other ODM methods provides an additional con- 
firmation of our conjecture about the sign of the discon- 
tinuity. 

VIII. CONCLUSIONS 

Let us briefly summarize the findings of the current nu- 
merical study of properties of the cubic anharmonic os- 
cillator. For g > 0, the Hamiltonian ([I} is VT symmetric 
and has a real positive spectrum. The energy eigenval- 
ues have divergent, Borel summable expansions in pow- 
ers of the coupling constant g. The imaginary parts 
of the eigenvalues on the cut on the real negative axis 
are positive. Starting from the VT symmetric case and 
developing a strong coupling expansion in the variable 
\J\ = .g 1 / 5 , one sees that eigenvalues are univalued func- 
tions in some neighbourhood of the origin in the £ plane, 
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TABLE VI: Values for £ [qc, (x), obtained by various ODM summation methods, on the real neg- 
ative axis, at transformation order 150. Method ()42a[) and the continued fraction approximation 
(|99[) provide the best numerical convergence. 



X 



_ 2 4/5 



4.741 



Method of Eq. (|42al) 
Method of Eq. (l42bl) 
Method of Eq. (f42c|) 
Continued fraction 



no convergence 
0.42 ±0.20 - (0.39 ±0.17) i 
0.38(1) - 0.36(8) i 
0.3898(5) - 0.3644(3) i 



0.19575 08157(1) 

0.195(8) 

0.195(5) 

0.19575 08157161719(6) 



X 



-5 4/5 = -0.275 . 



-(21.6) _4/5 = -0.085. 



0.28269 92581 93274 90989 9(0) 0.34215 80186 19340 42140 767(3) 

0.28269 93(0) 0.34215 80(2) 

0.2827(1) 0.34215(6) 

0.28269 92581 93274 90989 90(1) 0.34215 80186 19340 42140 767(6) 



Method of Eq. (|42a|) 
Method of Eq. (|42b| 
Method of Eq. ((42c)) 

Continued fraction 



-1,6 -f,4 -1,2 



-0,8 -0,6 -0,4 -0,2 



FIG. 3: Aoi = (£ [ gs - £ 1 es)7 4 as a function of X for 
— 1.85 < x < —0.85 [the value at the origin is Aoi(x = 0) = 
0.2263073... 1. 



FIG. 4: Soi = (£ [ gs ±£ , es)/2 as a function of X for -1.85 < 
X < —0.85 [the value at the origin amount to Sbi(x = 0) = 
0.8482634... 1. 



as a unique 



with a simple pole at z = = g l ' 5 
singularity. 

In Sec. IIIIl we have discussed various possible sum- 
mation methods for the determination of energy eigen- 
values of the cubic potential. Notably, we explore the 
numerical properties of three such methods here, given 
in Eqs. (|42a|) . (|42h>|) and (|42c[) . The latter was used in our 
previous paper [4J . Summing the perturbative expansion 
of the ground state energy E(g) by various implementa- 
tions of the ODM method, we uncover additional proper- 
ties. From the apparent convergence of the ODM method 
for a = 5/2 [see Eq. (|42c[) ]. we conclude that in the vari- 
able z the function E(g(z)) has no other singularity in 



the sector |argz| < 7r/4. In the variable x = 9 



-4/5 



energy (x) of the Hamiltonian ([5]) also is a real an- 
alytic function in a cut plane with a cut along the real 
negative axis. 

The ODM summation method (|42a[) converges very 
well in the strong coupling regime and thus allows us to 
determine a number of terms of the large g, thus small x 
expansion. The precision can be further improved by the 
convergence acceleration algorithm (|64[) and, finally, by a 
combined calculation of the ground state and first excited 
energies. The first 20 terms are given in Eq. ([92]) . A 
direct analysis of the strong coupling expansion indicates 
the existence of a square root singularity located at x = 
Xc ~ —1-351 [g w 0.687 exp(5i7r/4) in the second sheet]. 
These properties are confirmed by the direct calculation 



this translates into the property that the ground state of the difference (Eqs ~ Ees) , which vanishes linearly at 
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the point Xc = -1.3510415966(3), as given in Eq. (|9T1) . 

The strong coupling expansion, in turn, can be 
summed by an expansion in a continued fraction 
(Sec. IVIIB[) . The calculated coefficients of the contin- 
ued fraction display an asymptotic behavior consistent 
with a square root singularity. The extrapolated contin- 
ued fraction yields results in remarkable agreement with 
the more direct ODM calculations. 

As a function of z, E^ qc \xi z )) is invariant under the 
rotation z i— > zexp(i7r/2). It is thus entirely defined by 
its values for — 7r/4 < arg(z) < +7r/4, that is, — 57r/4 < 
ai 'g(<?) < +57r/4. As a consequence, as a function of z, 
E has cuts only on the lines argz = ti-k/A where level 
merging can occur. Note that, from the point of view of 
Eq. (|T7|) . for x < the Hamiltonian is still VT symmetric 
and, thus, the singularity at \ = Xc corresponds to a kind 
of spontaneous VT symmetry breaking. 

All coefficients of the associated expansion in a con- 
tinued fraction are found to be positive as long as we 
can estimate them with enough precision, that is, up to 
order 27 (corresponding to a strong coupling expansion 
up to order 28), and seems to converge toward a positive 
value consistent with the square root singularity of the 
first level merging. This gives a very strong argument 
that the eigenvalue of the strong coupling Hamiltonian 
Im£ , [ qc l(— \x\ + i0) is negative on the cut. This conjec- 
ture is also consistent with the large x behaviour, the 
Cauchy representation, the coefficient of the square root 
singularity. However, this property cannot be shared by 
all eigenvalues and may even be unique to the ground 
state. 

To go beyond this study, one would have to apply these 
ODM summation methods more directly to the spectral 
equation, in order to eliminate all level merging singular- 
ities (for the pure cubic potential see Ref. |2J|). Finally, 
we would like to emphasize the empirical evidence gath- 
ered during the current study, for the robustness of the 
ODM summation methods (|42a|) . (|42b[) and (|42c|) . which 
applied to different functions with varying implementa- 
tion always give consistent results. 



form, 

E(g)=gPJ2tn9- n/a , 9 -> oo , (Al) 

n=0 

a property shared by the example (|17p we discuss here. 
It is also shared by the quartic anharmonic oscillator [l3| 
and all x N perturbations to the quantum harmonic os- 
cillator. We then consider the conformal mapping (14T1) . 

This transformation maps the real positive g axis onto 
the finite A interval [0, 1]. For g — > oo, A — s- 1 and g has 
an expansion at A = 1 of the form 

<r 1/a = X>«(! - A ) n+1 ( A3 ) 

?i=0 

with A = p~ 1/a . The function 

0(A) = (1 - X) al3 E(g(X)) (A4) 
then has a Taylor series expansion at A = 0, 

0(A) = J2 ^ X n , (A5) 

n=0 

as well as at A = 1, 

0(A) = X>n(l-A) n . (A6) 

n=0 

with ipo = eq , where eo is the coefficient defined 
in (jAlj) . This last property explains, to a large extent, 
the good convergence of the method even for g — > oo. 



Appendix B: Large g expansion: a few terms 
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Appendix A: ODM and strong coupling expansion 

We assume that we know a perturbative expansion for 
E(g) in powers of the coupling constant g, and that E{g) 
also has a strong-coupling asymptotic expansion of the 



In order to determine the successive terms of the large 
g, small x expansion we do not have to differentiate nu- 
merically the values of E{g) calculated at g large but 
finite. In the ODM method, the various terms have ex- 
plicit analytic forms and we give here the first terms. 
First, we set 

E(g) = <f>(X)(l - X)~^ 

= 0(A) ^ A -^(j)%(A). (Bl) 
Then we write the relation (|A2[) as 

A 1 /". (B2) 
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Setting (p/g) 



l/a 



z, we expand s = 1 — A in powers of 



i, so that s — 1 — A = z 



z' + 



. The function "0(A) 



can then be expanded successively in powers of ,s and z: 
^(A) = 0(l)-^'(l) s + iV"(l)s 2 + --- 



= p f ) ( eo + p -V<* eiZ + p -V<* e2Z 2 



, (B3) 



where the e n coefficients are given in Eq. (|A1[) . The first 
terms yield 



e = 0(1) p-f, 



(B4) 



£3 



A shift of the variable x —> x — \ u then results in 

6 



- «V5 



Id 2 n u 2 u 3 

'- — ~ + X X H 

2 <9x 2 12 180 



,V5 



1 d 2 

2 <9x 2 



-4/5 



12' 



x+ im 9 



-6/5 



(C4) 



If one now writes the strong-coupling expansion of the 
iVth energy level in the form of Eq. (17) of Ref. [H, 



En(9) =5 1/5 Y, L ^9' 2K/5 



(C5) 



then the terms [see Eq. (|C4[) ] 



1 



(C6) 



Appendix C: Numerical evidence for strong coupling 



In Ref. [22[, we started from the Hamiltonian 

1 



2 dx 2 



1H2 + 2 X+ ^ X 



(CI) 



which entails both a change in the normalization of the 
coupling term and also a change in the complex phase of 



the coupling. The transformation x 
the Hamiltonian into the form 



9 



-1/10 



H 



9 



1/5 



We now define g = i 
H = g^ 



1 d 2 , 
-2^ + ^ 4 

" 5 / 2 and write 

i a 2 3 
2dx 2+ x ' 



2 a 5 / 2 



a; brings 



(C2) 



(C3) 



generate a convergent series in even powers of g 4 ^ 5 , and 
the term 



-6/5 



180 



(C7) 



in Eq. (|C4[) is the only term of odd power in g~ 2 ^ in 
Eq. (|C5p . From the point of view of the strong coupling 
expansion, this term Ljq.s = 1/108 represents a shift, 
independent of the quantum number N, of all levels of 
the cubic potential. This affords an explanation for the 
observations made in Table 2 of Ref. [22( , where on the 
basis of numerical evidence, it was conjectured that the 
terms of order K = 1 according Eq. (|C5[) vanish, and 
that the numerical value for the coefficient of the term of 
order K = 3 is 1/108 uniformly for all levels. 
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